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Abstract We investigate the connections between several recent methods for the discretization of aniso- 
tropic heterogeneous diffusion operators on general grids. We prove that the Mimetic Finite Difference 
scheme, the Hybrid Finite Volume scheme and the Mixed Finite Volume scheme are in fact identical up 
to some slight generalizations. As a consequence, some of the mathematical results obtained for each of 
the method (such as convergence properties or error estimates) may be extended to the unified common 
framework. We then focus on the relationships between this unified method and nonconforming Finite 
Element schemes or Mixed Finite Element schemes, obtaining as a by-product an explicit lifting operator 
close to the ones used in some theoretical studies of the Mimetic Finite Difference scheme. We also show 
that for isotropic operators, on particular meshes such as triangular meshes with acute angles, the unified 
method boils down to the well-known efficient two-point flux Finite Volume scheme. 

1 Introduction 

A benchmark was organized at the last FVCA 5 conference [HI in June 2008 to test the recently developed 
methods for the numerical solution of heterogeneous anisotropic problems. In this benchmark and in this 
paper, we consider the Poisson equation with homogeneous boundary condition 

-div(AVp) = / in fl, (1.1a) 
p = Q on9il, (1.1b) 

where il is a bounded open subset of R d (d> 1), A : — > M<j(R) is bounded measurable symmetric and 
uniformly elliptic (i.e. there exists £ > such that, for a.e. x £ Q and all £ £ R d , A(x)£ ■ £ > £l£| 2 ) an d 

feL 2 (n). 

The results of this benchmark, in particular those of [51 EE [HJ H2] seem to demonstrate that the behavior 
of three of the submitted methods, namely the Hybrid Finite Volume method [TH [T71 ITS] , the Mimetic 
Finite Difference method 0), and the Mixed Finite Volume method [TUl \H\ are quite similar in a 
number of cases (to keep notations light while retaining the legibility, in the following we call "Hybrid" , 
"Mimetic" and "Mixed" these respective methods). 

A straightforward common point between these methods is that they are written using a general partition 
of fl into polygonal open subsets and that they introduce unknowns which approximate the solution p 
and the fluxes of its gradient on the edges of the partition. However, a comparison of the methods is still 
lacking, probably because their mathematical analysis relies on different tool boxes. The mathematical 
analysis of the Mimetic method [2] is based on an error estimate technique (in the spirit of the mixed 
finite element methods). For the Hybrid and the Mixed methods [16, 10 . the convergence proofs rely 
on discrete functional analysis tools. The aim of this paper is to point out the common points of these 
three methods. To this purpose, we first gather, in Table [T] of Section [5J some definitions and notations 
associated with each method, and we present the three methods as they arc introduced in the literature; 
we also present a generalized or modified form of each of the methods. The three resulting methods are 
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then shown to be identical (Section [3]) and to inherit some of the mathematical properties of the initial 
methods (Section H}. Particular cases are then explored in Section [5j which show the thorough relations 
between this unified method and a few classical methods (nonconforming Finite Elements, Mixed Finite 
Elements and two-point flux Finite Volumes) and, as a by-product, provide a flux lifting operator for a 
particular choice of the stabilization parameters (but any kind of grid) . Finally, a few technical results 
are provided in an appendix. 

In this article, we only consider the case of a linear single equation , but the three methods we study 
have also been used on more complex problems, such as the incompressible Navier-Stokes equations 
[HI [TT] , fully non-linear equations of the p-Laplacian type [2 [T^] , non- linear coupled problems [7] , etc. 
However, the main ideas to apply these methods in such more complex situations stem from the study 
of their properties on the linear diffusion equation. The unifying framework which is proposed here has 
a larger field of application than 1|) ; it facilitates the transfer of the ideas and techniques used for one 
method to another and can also give some new leads for each method. 

2 The methods 

We first provide the definitions and notations associated with each method (in order for the readers who 
are familiar with one or the other theory to easily follow the subsequent analysis, we shall freely use one 
or the other notation (this of course yields some redundant notations). 





Mimetic notation 


Finite Volume notation 


Partition of in polygonal sets 




M 


Elements of the partition (grid cells) 


E 


K ( "control volume" ) 


Set of edges/faces of a grid cell 


dE, or numbered 
from 1 to kp] 


S K 


Edges/faces of grid cell 


e 


a 


Space of discrete p unknowns (piecewise approxi- 
mations of p on the partition) 


Q h 


Hm 


Approximation of the solution p on the grid cell 
K= K 


PE 


PK 


Discrete flux: approximation of J e — AVp ■ n E 
(with e = a an edge of E = K and n E = riK,a the 
unit normal to e outward E) 


F% 


Fk.o 


Space of discrete fluxes (approximations of 


X h 


T 



Table 1: Usual notations and definitions in Mimetic and Finite Volume frameworks. 



Remark 2.1 In the usual Finite Volume literature, the quantity Fk,o is usually rather an approximation 
of J a — AVp ■ Ta.K,a than j a — AVp • n^,,; we choose here the latter normalization in order to simplify 
the comparison. 

In all three methods, a natural condition is imposed on the gradient fluxes (called conservativity in Finite 
Volume methods and continuity condition in Mimetic methods): for any interior edge a (or e) between 
two polygons K and L (or E\ and £2), 

F K , a + F L ,a = 0, (or F Ei = -F% 2 ). (2.1) 

This condition is included in the definition of the discrete flux space T (or X h ). 
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The Mimetic, Hybrid and Mixed methods for all consist in seeking p G Hm (or Q h ) and F G T (or 
X h ), which approximate respectively the solution p and its gradient fluxes, writing equations on these 
unknowns which discretize the continuous equation Each method is in fact a family of schemes 

rather than a unique scheme: indeed, there exists some freedom on the choice of some of the parameters 
of the scheme (for instance in the stabilization terms which ensure the coercivity of the methods). 

2.1 The Mimetic method 

The standard Mimetic method first consists in defining, from the Stokes formula, a discrete divergence 
operator on the space of the discrete fluxes: for G G X h , VTV h G G Q h is defined by 

(VlV h G) E = ^Y,\e l \G%. (2.2) 
' ' i=i 

The space Q h of piecewise constant functions is endowed with the usual L 2 inner product 

\p,q\Q» = X! \ e \p e< i e 

and a local inner product is defined on the space of fluxes unknowns of each element E: 

k E 

[F e ,G e ]e = F E M E G E = £ M E . s , r F c gG^, (2.3) 

s,r=l 

where Me is a symmetric definite positive matrix of order k E . Each local inner product is assumed to 
satisfy the following discrete Stokes formula (called Condition (S2) in [2JE]): 

V£ G Q h , Vg affinc function. VG G X h : [(AVq) 1 , G] E + I <] {VlV h G) E dV = G E / qdT, (2.4) 

JE j=1 Je t 

where {{AVq) 1 )^ = J e A E 'Vq ■ n E dT, and A# is the value, assumed to be constant, of A on £ (if A is 
not constant on E, one can take equal to the mean value of A on £). 

Remark 2.2 Note that Condition of is needed in the convergence study of the method, but 

not in its definition; therefore it is only recalled in Section^ see (|4.49[) . 

The local inner products (|2.3[) allow us to construct a complete inner product [F, G] X h = J^Een, G]e, 
which in turn defines a discrete flux operator Q h : Q h — > X h as the adjoint operator of VTV h : for all 
F G X h and p G Q h , 

[F,g h p] xh = [p,VlV h F] Qh 

(notice that this definition of Q h p takes into account the homogeneous boundary condition (jl.lbjl ). Using 
these definitions and notations, the standard Mimetic method then reads: find (p, F) G Q h x X h such 
that 

VXV h F = f h , F = g h p (2.5) 

where //,, is the L 2 projection of / on Q h , or in the equivalent weak form: find (p,F) G Q h x X h such 
that 

VG G X h : [F, G] xh - \p, VlV h G] Qh = , (2.6) 
Vq G Q h : [DlV h F, q] Qh = [f hl q] Qh . (2.7) 
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The precise definition of the Mimetic method requires to choose the local matrices M.e defining the local 
inner products [•,•]£• It can be shown [3] (see also Lemma 16.21 in the appendix) that this matrix defines 
an inner product satisfying (|2.4p if and only if it can be written MbNb = M_e or equivalently 



where 



M.e is the kE x d matrix with rows (|ei|(a; ei — XE) T )i=i,k B i 
with x e the center of gravity of the edge e and Xe the center of gravity of the cell E, 

Ce is a Ue x (ks — d) matrix such that Im(C.E) = (Im(N.E)) , ^ 
with Ne the kE x d matrix with columns 
/ (A B )i-ng \ 
(N B )i= i for j = l,... ,d, 

(As)j being the j-th column of 



(2.9) 



(2.10) 



Ue is a (fc^; — d) x (kE — d) symmetric definite positive matrix. (2-11) 

Here, we consider a slightly more general version of the Mimetic method, replacing by a point Xe 
which may be chosen different from the center of gravity of E. We therefore take 

M E = -±-R E A E lR l + CeVeCI (2.12) 
\M 

where 

Re is the kE x d matrix with rows (|ej|(x ei — XE) T )i=i,k E i in -.on 

with x e the center of gravity of the edge e and xe any point in the cell E. 

The other matrices C E and U E remain given by (f2~TU|) and (j2~TT]) . The choice [ ([27T2j) . ((XT5]) ] of no 
longer gives, in general, an inner product [■, -]e which satisfies (|2.4j) . but it yields a generalization of this 
assumption; indeed, choosing a weight function we ■ E — > R such that 

/ w E (x) dx = \E\ and / xwe(x) dx = \E\xe, (2-14) 

we prove in the appendix (Section 16. ip that the matrix can be written (|2.12p with Re defined by 
(|2. 13|) if and only if the corresponding inner product [•, -]e satisfies 

VE Efl h , \/q affinc function, VG E X h : [(AVq) 1 , G] E + / q{VlV h G) E w E dV 

Je 

k E . (2.15) 

1 = 1 Je i 

Definition 2.3 (Generalized Mimetic method) The Generalized Mimetic scheme for (jl.lj) reads: 

Find (p,F) E Q h x X h which satisfies [(l2^. (f231) . p^ . (f2^ . (l2^ . <l2T3l^ ] . 

Its parameters are the family of points (xE)Een h (which are freely chosen inside each grid cell) and the 
family of stabilization matrices { < CE^E)Een h satisfying (|2.10jl and (|2.11|) . 
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2.2 The Hybrid method 

The standard Hybrid method is best defined using additional unknowns p a playing the role of approx- 
imations of p on the edges of the discretization of f2; if £ is the set of such edges, we let Hj^ be the 
extension of Hm consisting of vectors p = ((pk)k<em, (p<t)<t££)- It will also be useful to consider the 
space Hk of the restrictions pk = (pk, (Pa)<reS K ) to a control volume K and its edge of the elements 
p £ Hm . A discrete gradient inside K is defined for px £ Hk by 

' ' <t££k 

The definition of this discrete gradient stems from the fact that, for any vector e, any control volume K 
and any xk £ R d , we have 

|if|e= Yj W\e ■ {xa ~ x K )n K ,a (2.17) 

rr££ K 

where x a is the center of gravity of er and xk is any point of K . Hence the discrete gradient is consistent 
in the sense that if p& = ip{x^) and pk = iP(xk) for an affine function ^ on K, then V kPk = V?/> 
on K (note that this consistency is not sufficient to ensure strong convergence, and in fact, only weak 
convergence of the discrete gradient will be obtained). A stabilization needs to be added to the discrete 
gradient (|2.16p in order to build a discrete coercive bilinear form expected to approximate the bilinear 
form (it, v) t— > J n AVm • Vi> dx occurring in the weak formulation of (jl.ip . Choosing a point Xk for each 
control volume K (for instance the center of gravity, but this is not mandatory), and keeping in mind that 
Pk is expected to be an approximation of the solution p of at this point, a second-order consistency 
error term S k {pk) = (S K ,a{PK))ae£ K is defined by 



Sk,c{pk) =Pa- Pk - ^kPk ■ {x a ~ x K )- (2.18) 



Note that thanks to (|2"T7|) . 



J2 W\S K APK)nKa = 0, (2.19) 

and that Sk,<t{pk) — if p a — ip{x a ) and pk = tJj(xk) for an affine function i\) on K. 

The fluxes {Fk,o-)<ti=£ k on the boundary of a control volume K associated with some p £ Hj^ are then 

defined by imposing that 

Vit £ M , Vq K = (qK, (q*)ae£ K ) £ H K ■ 

V \a\F K AlK - q„) = \K\A k VkPk ■ Vjr?jr + V a K><7 -j^-S k ,<t(j>k )S K ,„{q K ) (2.20) 

ct£L£k (y^tK 

where Ak is the mean value of A on K, d,K,u is the distance between xk and the hyperplane containing 
a and otK.a > 0. Note that Fk.o is uniquely defined by (|2.20j) . since this equation is equivalent to 

\ct\F k ,o = \K\A k VkPk ■ y K q K + V a K ,* -^-Sk^^SkAQk) (2.21) 

o-e£K 

where qK satisfies qK — q a = \ and qK — q a ' = for a' ^ a. To take into account the boundary condition 
(|l.lbp . we denote by if .mo = {p £ H^i such that p a = if a C 90} and the Hybrid method can then be 
written: find p £ Hm,o such that, with Fk,g defined by (|2.20|) . 

Vq£H Mfi : E W\FkA<Ik - q a ) = q K f f. (2.22) 

KeM ae£ K K&M " K 
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In particular, taking qK = for all K, and q a = 1 for one a and for the others in (|2.22p yields that F 
satisfies (|2.ip ; once this conservativity property is imposed by requiring that F G F, we may eliminate 
the q„ from (|2.22[) and reduce the Hybrid method to: find (p, F) G Hm,o x F satisfying (|2.20[) and 



VKeM : V \a\F K , a = [ /. 



(2.23) 



This last equation is the flux balance, one of the key ingredients of the finite volume methods. 



Let us now introduce a generalization of the Hybrid method for the comparison with the other methods. 
As previously mentioned, the stabilization term Sk is added in (|2.20[) in order to obtain a coercive scheme 
(the sole discrete gradient (VkPk)k£M does not allow to control p itself); the important characteristic 
of Sk is that it yields a coercive bilinear form, so that we may in fact replace (|2.20JI by the more general 
equation 

VK G M. , Vq K G H K : 
J2 \v\FkAik - q a ) = \K\A k VkPk -y K q K + R K,*^SK,*(pK)S K ^(qK) 

= \K\k K VKPK-VKqK + S K {q K ) T V> K S K {p K ), (2.24) 

where B|£ a a , are the entries of a symmetric definite positive matrix Bj^ (the superscript H refers to the 
Hybrid method). 

Definition 2.4 (Generalized Hybrid method) A Generalized Hybrid scheme for (jl.ip reads: 

Find {p,F) G H Mfi x T which satisfies [ $rnty . $rn% . $FZty . $TZty ]. 

Its parameters are the family of points (xk)k<£M (which are freely chosen inside each grid cell) and the 
family (Mk)keM of symmetric definite positive matrices. 

Remark 2.5 Another presentation of the Generalized Hybrid method is possible from (|2.22[) and (|2.24D 

by eliminating the fluxes: 

Find p £ Hm,o such that 

Vgei/A4,o : J2 \ K \^KVKPK-V K qK+ Sk(Qk) T KSk(pk) = J2 9k [ /■ (2-25) 

KeM kem KeM J K 



2.3 The Mixed method 

As in the Hybrid method, we use the unknowns p in Hm for the Mixed method (that is to say approximate 
values of the solution to the equation inside the control volumes and on the edges), and fluxes unknowns 
in T . However, contrary to the Hybrid method, the discrete gradient is no longer defined from p, but 
rather from the fluxes, using the dual version of (|2.17|) . that is: 

\K\e= ^ \a\e- n K ,a(xa - x K )- (2.26) 

<tG£k 

For F G F, we denote Fk = {FK,a)ae£ K its restriction to the edges of the control volume K and Tk is 
the set of such restrictions; if Fk G Fk, then ~vk{Fk) is defined by 

\K \A k vk(F k ) =~Y1 W\ f kA** ~ x k) (2-27) 

cr££ K 
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where, again, Aa is the mean value of A on AT, x a is the center of gravity of a and xk a point chosen 
inside K. Recalling that -Frt.ct is an approximation of t^t J a — AVp- n^ :(T , Formula (|2.26p then shows that 
vk(Fk) is expected to play the role of an approximation of Vp in K. 

The Mixed method then consists in seeking (p,F) G Hm,o x T (recall that F G T satisfies (|2.ip . and 
we impose p a = if a C 9f2) which satisfies the following natural discrete relation between p and this 
discrete gradient, with a stabilization term involving the fluxes and a positive parameter v > 0: 

VAT eM,Vae£ K : Pa- Pk = v K (F K ) ■ (x a - x K ) - v&iam{K)F K , a (2.28) 
and the flux balance: 



VKeM : J2 \ a \ FK <- = 




(2.29) 



Multiplying, for any Ga 6 Tk, each equation of (|2.28p by \o\Gk,<j and summing on a G £a ; we obtain 
VA' G .M , VGa G ^if : 

]T (p A - -Pa)HGjr,a = |A-|vjf (JVf) ■ Aav a -(Ga) + ]T i/diam(A') MAa, ct Ga, ct . 

a££ K a££ K 

The term X)cre£ K ^diam(lf )|<7|i<A' ) a-GA',<7 in this equation can be considered as a strong stabilization term, 
in the sense that it vanishes (for all Gk) only if Fk vanishes. We modify here the Mixed method by 
replacing this strong stabilization by a weaker stabilization which, as in the Hybrid method, is expected 
to vanish on "appropriate fluxes". To achieve this, we use the quantity 

TkAFk) = F K ,a + A k v k (Fk) ■ n K>rr , (2.30) 

which vanishes if {Fk,o)o^£ k are the genuine fluxes of a given vector e. Then, taking to be a 
symmetric positive definite matrix, we endow Tk with the inner product 

{Fk,G k )k = \K\v K {F K )-K K v K {G K )+ rA 1*,„>TkAFk)T k ^{G k ) 

a,a'££ K 

= \K\v K (F K ) ■ AaVa(Ga) + T K (G K fM^T K (F K ) (2.31) 
and the stabilized formula (|2.28p linking p and F is replaced by 

VA' G M , VGa G T k : (F K , G K ) K = ^ (p K - p a )\(r\G Kj(T . (2.32) 

a££ K 

We can get back a formulation more along the lines of (|2.28p if we fix a G £k and take Gk(o~) G Tk 
defined by Ga(c)o- = 1 and Ga(c)<t' = if a' ^ a: (|2.32p then gives 

Pa Pk = ~\K\v K {F K ) ■ A K VK(G K (a)) - ^T K {G K {o-)) T M^T K {F K ). 
I CT I \ a \ 

But |A'|AaVa(Ga(c)) = — IcrKxo- — xk) and thus 

Pa -PK = va(Aa) ■ (xa - x K ) - ^T K {G K (<j)fK T K {F K ) ■ 

This equation is precisely the natural discrete relation (|2.28p between p and its gradient, in which the 
strong stabilization involving Fk,<? has been replaced by a "weak" stabilization using Tk(Fk)- 

Definition 2.6 (Modified Mixed method) A Modified Mixed scheme for (ji.ip reads: 

Find (p,F) G H M ,o x T which satisfies [ $Z^ , $Z^ , $^ , $]n$ , $]E$ \. 

Its parameters are the family of points (xk)k£M (which are freely chosen inside each grid cell) and the 
family (R k i )k£M of symmetric definite positive matrices. 
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3 Algebraic correspondence between the three methods 



We now focus on the main result of this paper, which is the following. 

Theorem 3.1 (Equivalence of the methods) The Generalized Mimetic, Generalized Hybrid and Mo- 
dified Mixed methods (see DeHnitions \2.Sl \2.4\ and \2.b}) are identical, in the sense that for any choice of 
parameters for one of these methods, there exists a choice of parameters for the other two methods which 
leads to the same scheme. 

The proof of this result is given in the remainder of this section, and decomposed as follows: for comparison 
purposes, the Generalized Mimetic method is first written under a hybridized form in Section 13.11 then, 
the correspondence between the Generalized Mimetic and the Modified Mixed methods is studied in 
Section 13. 2\ finally, the correspondence between the Generalized Mimetic and the Generalized Hybrid 
methods is carried out in Section [3.31 



3.1 Hybridization of the Generalized Mimetic method 

Although the edge unknowns introduced to define the Generalized Hybrid and Modified Mixed methods 
may be eliminated, they are in fact essential to these methods; indeed, both methods can be hybridized 
and reduced to a system with unknowns {pa)ae£ on ly- I n order to compare the methods, we therefore 
also introduce such edge unknowns in the Generalized Mimetic method; this is the aim of this section. 
Let E be a grid cell and e £ dE be an edge. If e is an interior edge, we denote by E the cell on the 
other side of e and define G(e,E) £ X h by: G{e,E)% = 1, G(e,E) e ~ = -1 and G{e,E)%, = if e' ^ e. 

We notice that G(e, E) = —G(e, E) and, using G(e, E) in (|2.6j) . the definitions of VIV h and of the inner 
products on X h and Q h give \e\(p E p g ) = [F,G{e,E)] E + [F,G(e,E)} 5 = [F,G(e,E)] E - [F,G(e,E)] s . 
It is therefore natural to define p e (only depending on e and not on the grid cell E such that e C dE) by 

V£ £ fl h , Ve £ dE : \e\(p E - p e ) = [F,G(e,E)] E . (3.33) 

This definition can also be applied for boundary edges e C <90, in which case it gives p e = (thanks to 

(ED). 

We thus extend p £ Q h into an element p £ Hj^ q having values inside the cells and on the edges of the 
mesh. Denoting, as in Section |2"UI T E the space of restrictions G E to the edges of E of elements G £ X h 
and writing any G E £ T E as a linear combination of (G(E, e)) eS aE, it is easy to sec from (|3.33[) that the 
Generalized Mimetic method [ (|2.6p . (|2.7[) ] is equivalent to: find (p,F) £ -Ha^o x X h such that 

V£ £ n h , \/G E £ T E : [F e ,Ge]e = \ e \(P E -Pe)G E (3-34) 

e£dE 

and 

VE £ Q h : ]T \e\F% = f f. (3.35) 



3.2 Proof of the correspondence Generalized Mimetic <-> Modified Mixed 

The simplest comparison is probably to be found between the Modified Mixed and Generalized Mimetic 
methods. Indeed, we see from [ ()2.29|) . ()2.32[) ] and [ (|3.34|) . (|3.35j) ] that both methods are identical provided 
that, for any grid cell K = E, the local inner products defined by (|2.31|) and (|2.12|) are equal: (-,-)k = 
[•, -] E ,. We therefore have to study these local inner products and understand whether they can be 
identical (recall that there is some latitude in the choice of both inner products). 
Let us start with the term \K\vk(Fk) • ^-k^k{Gk) in the definition of (•, -)k- Thanks to (|2.27p . 

\K\v K (F K ) ■ A K v K (G K ) = tL[J2 A k(I ct K^ -*k))Fk,<, ) ■ ( ^ W\(^ - x K )G K ,a) ■ 
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But, with the definition (|2.13[) of Re, ^2 ae£K \(?\(xcr — xk)Fk,o = R^j-Fir and thus 

\K\v K (F K ) ■ A K v K {G K ) = -L (A^rIFk) • (R^Gjc) = G£ (^R^A^R^ F K . 

Hence, the term |if |v^(-) • Akvk(-) in the definition of (•, -)k is identical to the first term p^-R^A^R^ 
in the definition of the matrix Me of [-, -]e (sec (|2.12[) ). Therefore, in order to prove that (■, -)k — [•, -]e, it 
only remains to prove that, for appropriate choices of Ce, Ue and B^, we have for any (Fk, Gk) & 3~ K : 

T k {G k ) T KTk{F k ) = GlC E V E ClF K (3.36) 

(see (pOTj) and ([2~12|) ): in fact, this is the consequence of Lemma 16 . 31 in the appendix and of the following 
lemma. 

Lemma 3.2 The mappings T K : K kE — > K kE and : H kE — > R fe B- rf have the same kernel. 

Proof of Lemma [3721 

We first prove that: 

i) Iiii(Ne) C Ker(Tft-) i.e. Tjc,cr((^E)j) = for all a 6 £k arid all j = 1, . . . , d (which also amounts 
to the fact that the lines of Tk are orthogonal to the vectors (N E )j). 

ii) dim(Im(TR-)) = Ue — d, and therefore dim(Ker(Tj<-)) = d. 

Item i) follows from (|2.27j) and (|2.26|) : we have A^v/f ((Ne)j) = — r^r J2ae£ K \ a \(^- K )o ' n K,a{x a — xk) = 
-(A K )j and thus T K ,a((N E )j) = (A K )j ■ n K , a - (A K )j • n A ',<r = 0. 

In order to obtain Item ii), we first remark that He — d is an upper bound of the rank of Tk since the 
lines of Tk arc in the orthogonal space of the independent vectors ((N E )j)j=i....,d Q)- Moreover, (|2.27[) 
shows that the rank of v A : H, kE — > R, kE is the rank of the family (x a — xk)<t<e£ k , that is to say d, and 
the kernel of v A therefore has dimension &e — d; since Tk = Id on this kernel, we conclude that the rank 
of Tk is at least k E — d 7 which proves ii). 

These properties show that Ker(Tft-) = Im(Ns) = (Ihi(Ce)) = ker(C^), and the proof is complete. ■ 

The comparison between Tk(Gk) T ^kTk(Fk) and G^CeUeC^Fk- is now straightforward. Indeed, 
let (Ce,Ue) be any pair chosen to construct the Generalized Mimetic method; applying Lemma 16. 3( 
thanks to Lemma 13.21 with A = C|J, B = Tk and {-,-}nk B -d the inner product on R fe E-a corre- 
sponding to Ue, we obtain an inner product {-,-};r*=.e on R feE such that {Tk(Fk),Tk(Gk)}-rj°e = 
{C^Fk, C^Gif} R /o B -d = G k CeV e'C'e-Fk', this exactly means that, if we define B^ as the matrix of 
{") '}R fc -B ; P-36p holds. Similarly, inverting the role of and Tk when applying Lcmma HTIfl for any Bf/ 
used in the Modified Mixed method we can find Ue satisfying (|3.36|) and the proof of the correspondence 
between the Generalized Mimetic and Modified Mixed methods is concluded. 



3.3 Proof of the correspondence Generalized Mimetic ^ Generalized Hybrid 

To compare the Generalized Mimetic and Generalized Hybrid methods, we use a result of [H [23] which 
states that the inverse of the matrix Me in (|2.12[) can be written 

M £ ' = We = ^NeA^nI + DeUeD| (3.37) 

where Be is a kE X (&e — d) matrix such that 

Im(D>E) = (Im^E))^ (3.38) 

6 Let us notice that the independence of ((Ni;)j)j = i can be deduced from the independence of the columns of A^: 
thanks to 1 12. 261 1 , any non-trivial linear combination of the (Ng ) j gives a non-trivial combination of the columns of A/f . 
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and Ve is a symmetric definite positive (k E — d) x (k E — d) matrix (note that the proof in [4ll23| assumes 
xe to be the center of gravity of E, i.e. M. E to satisfy (|2.8[) . but that it remains valid for any choice 
of xe, i-e. for any matrix Mb satisfying (|2.12[1 ). This result is to be understood in the following sense: 
for any (C E ,U E ) used to construct M E by (|2~T2")) . there exists (B E ,V E ) such that W E defined by (j3~3T)) 
satisfies Wb = , and vice-versa. 

For pe — (p E , {Pe)e£dE) G He, we denote by Ve the vector in R fcjs with components (\e\(pE — Pe))eedE- 
The relation (|3.34j) can be re-written M. E F E = V E , that is to say F E = Wb?b, which is also equivalent, 
taking the inner product in R fcfi with an arbitrary Q E (built from a q E G H E ), to 

-iE G n h ,Vq E eH E : J2 Mfe - Q*) f e = Qe^eVe (3.39) 

eedE 

The Generalized Mimetic method [ (|3.34p . (|3.35|) ] is therefore identical to the Generalized Hybrid method 
[ (f2T23]) . (f2T24|) ] provided that, for all E = K G fl h , 

V(p B) £3) G H E : QIWbT'b = \K\A K V K p K ■ V K q K + S K {q K ) T M^S K {p K ). (3.40) 

As in the comparison between the Generalized Mimetic and Modified Mixed methods, the proof of (|3.40|) 
is obtained from the separate study of each term of (|3.37D . 

First, by the definition (|2.10p of the matrix Ne and the definition (|2.16[) of VkPk, we have (NU'P.e)^ = 
Y^eedE l e l( A s)j • n|(p B - p e ) = -\K\(A K )j ■ VkPk for all j = 1, . . . , d, that is to say, by symmetry of 
A, NlV E = ~\K\A K V K p K . Hence, 

Q| ^N^A^N^ V E = \K\{k K V K q K ) T K K \k K V K p K ) = \K\K K V K p K ■ V K q K . 

The first part of the right-hand side in (|3.40j) thus corresponds to the first part of in (|3.37|) . 
and it remains to prove that (with appropriate choices of Be, U e and Of?-), for all (pe,Qe) £ He, 
QeS> e IS e B%V e = S K {qi<) T ^K Sk(pk)- Let us notice that, defining L K : R kE R kE by 

L K (V) = (y-fVa - D K V ■ (x a - x K )) with D K V = 4t V V a n K .„, (3.41) 

this boils down (letting V = Ve and V = Qe) to proving that 

V(V,V')eH kE : (V') T D E V E DlV = L K (V') T M K L K {V). (3.42) 

As previously, this will be a consequence of Lemma T6. 31 and of the following result. 

Lemma 3.3 The mappings Lk ■ R fcfi — > R fcfi and : R feE — > R fe E~ d fta^e i/ie same kernel. 

Proof of Lemma 13.31 From (|3.38|) . we get that Ker(D|J) = Im(R£). Hence it remains to prove that 
lm(R E ) = Ker(L K ) 

Let us first prove that Im(R£) C Kei(L E ). The j'-th column of M. E is (M. E )j = (|cr|(a;^. — x k:))^gS k (* nc 
superscript j denotes the j-th coordinate of points in R d ). Thus, for any e G R d , by (|2.26p . 

D K • e = — M e ■ n K,a(xi - x 3 K ) = e J , 

1 1 c^Sk 

which means that is the j-th vector of the canonical basis of R d . We then have DK(^ E )j ■ 

{x a — Xk) = — %k an< ^ ^ nUS 

{L K {(R E )j)) a = 4 - x 3 K - D K (R E )j ■ {x* ~ x K ) = 0; 
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this proves that the columns of R# are in the kernel of Lk, that is Iul(Re) C Ker(L#-). 
Next we notice that the rank of the mapping Dk '■ V £ R kE i— » DkV G R d (i.e. the rank of the fam- 
ily (T3.K,a)aeE K ) ls d an d that the mapping Lk is one-to-one on Ker(Dif). Hence dim(Im(Ljf)) > 
dim(Ker(Z?#-)) = &e — d, and therefore dim(Ker(Lx)) < d. But since Im(Re) C Ker(L^) and 
dim(Im(R£;) = d (the rank of the rows of R_e), we thus get that Im(RE) = Ker(Xx)- ■ 

From Lemmas 16.31 and 1331 we deduce as in Section l3~21 that for any choice of (D_e,Ue) there corresponds 
a choice of BjJ such that (|3.42|) holds, and vice-versa, which concludes the proof that the Generalized 
Mimetic method is identical to the Generalized Hybrid method. 

Remark 3.4 These proofs and Remark \6.4\ show that one can explicitly compute the parameters of one 
method which gives back the parameters of another method. Notice also that the algebraic computations 
required in Remark \6.4\ to obtain these parameters are made in spaces with small dimensions; the cost of 
their practical computations is therefore very low. However, to implement the methods, it is not necessary 
to understand which (Ce,Ue) or (De,Ue) corresponds to which B|/ orM 1 ^, since the only useful quan- 
tities are CeUeC^, OeUeB^, Tft-(-) T B^T/f (•) and Lk{-) t 1$>kLk{-), and the relations between these 
quantities are trivial (see (pHB]) and (f3~15|) ). 



4 Convergence and error estimates 

We showed in Section [3] that the three families of schemes which we called Generalized Mimetic, Gener- 
alized Hybrid FV and Modified Mixed FV are in fact one family of schemes, which we call the HMMF 
(Hybrid Mimetic Mixed Family) for short in the remainder of the paper. We now give convergence and 
error estimate results for the HMMF method. 



4.1 Convergence with no regularity 

In this section, we are interested in convergence results which hold without any other regularity assump- 
tion on the data than those stated in Section [TJ We therefore consider that A is only bounded and 
uniformly elliptic (not necessarily Lipschitz continuous or even piecewise continuous), that / € L 2 (fl) 
and that the solution to only belongs to H^fl) (and not necessarily to H 2 (fl)). 
We study how existing results, previously established for the Hybrid scheme, can be extended to the 
HMMF. In [16], we proved the L 2 convergence of the "standard" Hybrid method for a family of partitions 
of such that any cell K is star-shaped with respect to a point Xk and such that there exists 9 > 
satisfying, for any partition of the family, 



da a diam(in \ 
. ' . max 



max — — , max < 6, (4.43) 

<7££i«. d L . a KeM d K ,a / 

where £i n t denotes the set of internal edges of the mesh and M. a the set of cells to which a is an edge. We 
now show how the convergence of the HMMF may be deduced from an easy extension of (TBI Theorem 
4.1, Lemma 4.4] provided that: 

• in the Generalized Mimetic formulation [ ([^ . ((T3T) . (r^ . (r^ . (r^T^ . ((^l31) ]. 

there exist s* > and S* > 0, independent of the mesh, such that, 

for all cell K and all V = (V a ) a( z £K , . 

J2 W\ d K,AVa) 2 <V T M K V <S* W\dK,a(V a ) 2 , 
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in the Generalized Hybrid formulation [ (|2~TB1) . ([2~TB1) , ((235)) . (|2~2^|) ]. using the notation (|3"H|) . 

(4.45) 



there exist s„ > and 5* > 0, independent of the mesh, such that, 
for all cell K and all V — {V a ) a £S K : 

**E ffe£je P-(Lk(V))1 < L K {V) T M K L K {V) < S.^e* P~( L k(V)) 



• in the Modified Mixed formulation [ (j2^ . (|23^ . (|2^ . (|23T|) . (|232l) ]. 

there exist s* > and 5* > 0, independent of the mesh, such that, 

for all cell K and all V = (K-) CT e£ K) } (4.46) 

**E CT6fe WkATkAv)) 2 < T K (v) T n%T K (v) < s*£ CT£fK W\d K AT K Av)Y 

The three conditions (|4.44p . ()4.45p and (|4.46[) are in fact equivalent (this is stated in the next theorem), 
and one only has to check the condition corresponding to the chosen framework. 

Theorem 4.1 (Convergence of the HMMF method) Assume that A : fl — > M<j(R) is bounded 
measurable symmetric and uniformly elliptic, that f 6 L 2 (fl) and that the solution to belongs 
to Hq(£1). Let 6 > be given. Consider a family of polygonal meshes of fl such that any cell K is 
star-shaped with respect to a point xk, oxid satisfying (|4.43|) . Then the three conditions (|4.44|) . (|4.45l) 
and (|4.46|) are equivalent. Moreover, if for any mesh of the family we choose a HMMF scheme such 
that one of the conditions (|4.44p , (|4.45[) or (|4.46p holds, then, the family of corresponding approximate 
solutions converges in L 2 (Q) to the unique solution of (|l.ip as the mesh size tends to 0. 

Proof of Theorem 14.11 

Let us first prove the equivalence between (|4.44p and (|4.45p . assuming (|4.44p to begin with. Using (|3.42p , 
we get that, for all V, 

V t B k U k ^kV = L K (V) T M K L K (V). 

Let us apply the above relation to Lk(V) defined by LK(V) a = \a\LK(V) a . From (|2.26p (see also 
(|2.17p ) and recalling the operator Dk defined in (|3.4ip . we easily get that DkLk(V) = 0, which provides 
Lk(Lk(V)) = Lk(V). Hence we get 

(L K {v)) T n K v K n T K L K {v) = l k (v) t Kl k {v). 

We then remark that N k Lk(V) = (once again from (|2.26p ) and therefore, using (|3.37p . 

L K (V) T M K L K (V) = {L K (V)fW K L K (V). (4.47) 



Let Ik be the diagonal matrix Ik = d\a.g((^J\a\dK.a)a££ K ))• Condition (|4.44j) applied to Gk = IrFk 
gives, for all Gk 6 Tk, 

s*G T K G K < GkI^MkIAGk < S*G K G K . 
This shows that the eigenvalues of I^M^I^ 1 are in [s*, £»], and thus that the eigenvalues of Ir-M^- Ir- = 
I k W k Ik belong to Translating this into bounds on {I K 1 L K {V)) T I K ^ kIk^^r^)), we 

deduce 

and (|4.45p follows from (|4.47p (with s* = -g- and S* = j-). Reciprocally, from (|4.45p . following the proof 
of Lemma 4.4] and setting V a = \(t\(pk —p<j) in that proof, we get the existence of c\ > and c 2 > 0, 
independent of the mesh, such that 



Cl E rry — v ° ^ v t ^kV < c 2 ]T tt^ v °- ( 4 - 
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Using Ik as before, we then get (|4.44[) with s* = ~ and S* = 

Let us turn to the proof of the equivalence between (|4.44[) and (|4.46|) . beginning by assuming that (|4.44p 
holds. Using m> , one has T K ({A K v K (F K )-nK,a)aee x ) = 0, and thus T K {T K (F K )) = T K {F K ); hence, 
noting that M. k Tk(Fk) = (once again thanks to (|2.26p ) and remembering (|2.12p . (|3 . 36[) . we see that 
(|4.44p applied to V = Tk(Fk) directly gives (|4.46p . The reciprocal property follows, in a similar way as 
(I4.48p . from a simple adaptation of classical Mixed FV manipulations (used for example in the proof of 
a priori estimates on the approximate solution), see |10j . 

Note also that obtaining (|4.44[) from (|4.45p or (|4.46p is very similar to [4j Theorem 3.6]. 

We can now conclude the proof of the theorem, that is to say establish the convergence of the approximate 
solutions: using (|4.45p . this convergence is a direct consequence of [16, Theorem 4.1] with a straightfor- 
ward adaptation of the proof of [H)J Lemma 4.4] . The convergence could also be obtained, using (|4.46p , 
by an easy adaptation of the techniques of proof used in the standard Mixed setting (see [ID])- ■ 

Mimetic schemes are usually studied under a condition on the local inner products usually called (SI) 
and which reads [2]: there exist s» > 0, S* > independent of the mesh in the chosen family such that 

k E k E 

VE G Q h , VG G X h : s.^\E\(G e i) 2 <[G,G] E < S^\ e K G e) 2 - ( 4 - 49 ) 

i=l i=l 

The mesh regularity assumptions in [21 H] also entail the existence of C\ , independent of the mesh, such 
that 

qXjdiam^- 1 < H Va G £ K , . . 

qjjdiam(A') < d K ,„, Vcr G £ K . V ' ' 

Under such mesh assumptions and since |u| < 2 d ~ 1 diam(L) whenever a G £l, it is easy to see that there 
exists 9 only depending on CJjjsuch that (|4.43p holds; still using (|4.50p . we see that the quantities |cr|dft- iCr 
are of the same order as |.f£T| and thus that (|4.49p implies (|4.44p (with possibly different s* and 5*). We 
can therefore apply Theorem 14. II to deduce the following convergence result under the usual assumptions 
in the mimetic literature (except for the regularity assumptions on the data). 

Corollary 4.2 (Convergence under the usual mimetic assumptions) We assume that A : O — > 
Md(R-) is bounded measurable symmetric and uniformly elliptic, that f G L 2 (fl) and that the solution 
to p.ip is in Hq(Q,). Consider a family of polygonal meshes of fl such that any cell E is star-shaped 
with respect to a point xe, and such that (|4.50p holds (with Chn independent of the mesh in the family). 
We choose local inner products satisfying Condition (SI) of ^ (i.e. (|4.49p with s*, S* not depending on 
the mesh) and (|2.15p . Then, the family of approximate solutions given by the corresponding Generalized 
Mimetic schemes converges in L 2 (fl) to the solution of (jl.ip as the mesh size tends to 0. 
In particular, taking xe as the center of gravity xe of E, (|2.15p reduces to Condition (S2) of JMj (that 
is to say (|2.4p ), and the family of approximate solutions given by the corresponding "standard" Mimetic 
schemes converges in L 2 (fl) to the solution of (JTTTJ) as the mesh size tends to 0. 

Remark 4.3 (Compactness and convergence of a gradient) The proofs of the above convergence 
results rely on the use of the Kolmogorov compactness theorem on the family of approximate solutions, see 
\10[ Lemma 3.3] or \lb\ Section 5.2]. In fact, this latter study shows that for p G [l,+oo), any family of 
piecewise constant functions that is bounded in an adequate discrete mesh- dependent W 1,p norm converges 
in L P (VL) (up to a subsequence) to a function of W^' p (f2). The regularity of the limit is shown thanks to 
the weak convergence of a discrete gradient to the gradient of the limit \lb\ Proof of Lemma 5.7, Step 1]. 
Note that this compactness result and the construction of this weakly converging gradient are completely 
independent of the numerical scheme, which is only used to obtain the discrete W 1,p estimate and the 
fact that the limit is indeed the solution of (|l.ip . 

Moreover, it is possible, from the approximate solutions given by the HMMF schemes, to reconstruct 
a gradient (in a similar way as 1761 (22)-(26)]) which strongly converges in L 2 (fl) to the gradient of 
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the solution of (jl.ip , see 1761 Theorem 4..1]. One can also directly prove that the gradient Vk> defined 
from the fluxes by (|2.27p . is already a strongly convergent gradient in L 2 (fl), see in fact, the strong 
convergence of this gradient is valid in a more general framework than the one of the methods presented 
here, see ]T$. 

4.2 Order 1 error estimate 

Let us first consider the original Mimetic Finite Difference method (|2.2[) (|2. 1 1() . This method is consistent 
in the sense that it satisfies the condition (S2) of [2J. Under the assumptions that: 

• Condition (Ml) of [2] on the domain £1 (namely Q is polyhedral and its boundary is Lipschitz 
continuous) and Conditions (M2)-(M6) of (corresponding to (M1)-(M4) in g]) hold; 

• The stability condition (SI) of [SJ d] holds; this condition concerns the eigenvalues of the matrix 
Ue in (|2.8[) . and are connected with the discretization, 

• A € W hoc (n) dxd and p £ H 2 (Q), 

then Theorem 5.2 in [5J gives an order 1 error estimate on the fluxes, for an adequate discrete norm; 
moreover, if fl is convex and the right-hand side belongs to iJ 1 (fi), an order 1 error estimate on p in the 
L 2 norm is also established. 

By the equivalence theorem (Theorem 13. ip . this result yields similar error estimates for the Generalized 
Hybrid and Modified Mixed methods in the case where the point xk is taken as the center of gravity xk 
of K. 

In the case of the original Hybrid method (|2.22p , an order 1 error estimate is proved in [TB] only under the 
assumption of an homogeneous isotropic medium, that is A = Id, and in the case where the solution to 
(|l.ip also belongs to C 2 (f2) (but with no convexity assumption on the domain f2). As stated in Remark 

4.2 of [16], the proof is readily extended if the solution is piecewise H 2 (in fact, in [16] the situation is 
more complex because the Hybrid scheme is considered in the more general setting of the SUSHI scheme, 
which involves the elimination of some or all edge unknowns; see Section I5"H1 below) . 

4.3 Order 2 error estimate 

Under an additional assumption of existence of a specific lifting operator, compatible with the considered 
Mimetic Finite Difference method, a theoretical L 2 -error estimate of order 2 on p is proved in [2J. A 
condition on the matrix which ensures the existence of such a lifting operator is given in [3]: in 
particular, if the smallest eigenvalue of Ve is large enough (with respect to the inverse of the smallest 
eigenvalue of C^Cs and to the largest eigenvalue of a local inner product involving a generic lifting 
operator), then the existence of a lifting operator compatible with Me, and thus the super-convergence 
of the associated Mimetic scheme, can be proved. 

Regarding the consequence on the HMMF method, this means that, if xk is the center of gravity of K 
and if the symmetric positive definite matrices (U£;)_Bgo h or (R k )k£M or (Mk)kgm ar e "large enough", 
then the approximate solution p given by the corresponding HMMF method converges in L 2 with order 
2. It does not seem easy to give a practical lower bound on these stabilization terms which ensure that 
they are indeed "large enough" for the theoretical proof; however, several numerical tests (both with 
Xk the center of gravity of K, or xk elsewhere inside K) suggest that the HMMF methods enjoy a 
superconvergence property for a wider range of parameters than those satisfying the above theoretical 
assumptions. 

5 Links with other methods 

We now show that, under one of its three forms and according to the choice of its parameters, a scheme 
from the HMMF (Hybrid Mimetic Mixed Family) may be interpreted as a nonconforming Finite Element 
scheme, as a mixed Finite Element scheme, or as the classical two-point flux finite volume method. 
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5.1 A nonconforming Finite Element method 



In this section, we aim at identifying a hybrid FV scheme of the HMMF with a nonconforming Finite 
Element method. Hence wc use the notations of Definition 12.41 and Remark 12.51 For jury K G M and 
o G £k, we denote by A^ jCT the cone with vertex xk and basis a. For any given p G Hm we define the 
piecewise linear function p : Q — > R by: 



VA" G M , Vcr G £ K , for a.e. x G A K ,c ■ p( x ) = Pk + (vkPk + ^&,<rfe)njf,aj ■ (x - x K ) 

where (3k > (note that, since tik & • (x<7 — x/f) = dx,a, in the particular case where (3k = 1 we have 
p(x a ) = p a ). We let Hm = {p , p£ Hm} and we define the "broken gradient" of p G Hm by 

Vif G , Vcr G £a- , for a.e. a; G Ajr i<T : Vp{x) = Vp(x) = V ' kPk + -y^-S K .a(PK)n K .a- 

UK,cr 

We consider the following nonconforming finite element problem: find p G such that 

VqEH M ■ I A{x)Vp(x) ■ Vq{x)dx = [ f(x)q(x)dx, (5.51) 



where A(x) is equal to A^- for a.e. x G K. Then the above method leads to a matrix which belongs to the 
family of matrices corresponding to the Generalized Hybrid method (and thus also to the Generalized 
Mimetic and Modified Mixed methods). Therefore, since | A# i<7 | = ^ d J c '" , we get from (|2.19p : 

J l(x)Vp(x) ■ Vq(x) dx = \K\A K V K p K -V K q K + Yj S K (qK) T KSK(pK), 

TS /- K A 71 ' r~ \ A 



KeM KeM 

with 

Ik, a 

and, for a ^ a', 



kl(fe) 2 A 



^K,a,a' 



0. 



Note that this definition for fulfills (|4.45p under the regularity hypothesis (|4.43p . hence ensuring 
convergence properties which can easily be extended to this nonconforming Finite Element method (|5.5ip 
even though it is not completely identical to a scheme of the HMMF, since the right-hand sides do not 
coincide: in general, 



/ f(x)q(x) dx ^ Y Ik 

JQ J K 



KeM 

Indeed, this does not prevent the study of convergence since the difference between the two right-hand- 
sides is of order h. 



5.2 Relation with Mixed Finite Element methods 

In this section, we aim at identifying a particular scheme of the HMMF, under its mixed version with a 
Mixed Finite Element method. Let us first recall the remark provided in [21 Section 5.1] and Example 
1]: on a simplicial mesh (triangular if d = 2, tctrahcdral if d — 3), the Raviart-Thomas RTO Mixed Finite 
Element method fulfills properties (|4.49p and (|2.4p . and it is therefore possible to include this Mixed 
Finite Element scheme in the framework of the HMMF. Our purpose is to show that our framework also 
provides a Mixed Finite Element method on a general mesh (note that, even in the case of simplices, this 
method differs from the Raviart-Thomas RTO method). We use here the notations provided by Definition 
12.61 Still denoting A^-.er the cone with vertex xk and basis cr, we define, for F G T , 

VA' G A4 , Vcr G £k 7 Vx G A K , ff : F K , a (x) := -A K v K (F K ) + T K , a (F K )^^. (5.52) 

UK, a 
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If x belongs to the interface d A#- j(T ndAx^a-' between two cones of a same control volume, we have 
(x — xk) • H(9Ak ,naA K , = and thus the normal fluxes of Fk,<t are conservative through such interfaces; 
moreover, for all a G 8k and all x G a, we have (a;-%)'% i( , = <i.R- ;(T and thus, by (|2.30[) . Fk,<t{x)-t\k,<t = 
Fk,o\ since the elements of T satisfy (|2.ip . these observations show that the function F, defined by 

VK G M , Ver e Ek , for a.e. x G A KjJ : F(x) = F K ^(x), (5.53) 

satisfies F G H^wi^l). Noting that 

S 1 ^ = £ W^)^^ - (5.54) 

(thanks to (12.27p . (|2.30p and (|2.26p ). we have, with A again denoting the piecewise-constant function 
equal to in K, 

f A(x)~ 1 F(x) ■ G(x) dx = \K\v K {F K ) ■ A K v K (G K ) + V lK ..T K . a {F K )T K . a (G K ) (5.55) 



with 



iK.a = / A(x) — • — d.T > 0. 

JA k _ "A> d K ,a 



The right-hand side of (|5.55| defines an inner product {-,-)k which enters the framework defined by 
([2"UTj) , setting 

= lK,a, and M A J a a , = for a ± a', (5.56) 

which fulfills (|4.46p under the regularity hypothesis ()4.43p (hence ensuring convergence properties). There- 
fore the form [ (|2.6|) . ()2.7|) ] of the resulting HMMF scheme resumes to the following Mixed Finite Element 
formulation: find (p, F) G Hm x T such that 

VG G T : f A(x)- 1 F(x) • G(x) dx - j p{x)dWG{x) dx = , (5.57) 

Vq G H M : { q{x)divF{x) dx = [ q(x)f(x) dx , (5.58) 
Jn Jn 

where T = {F , F G J-}. Indeed, since | 1\k,o \ = and Y^ae£ K \&\iik,o- = 0, we have, for all F £ J 7 , 

j divF{x)dx= W\TkAFk)= Y, ^ 5 - 59 ) 

and (|5.58p written on the canonical basis of Hj^ is exactly (|2.29p . Summing (|2.32p on K, the terms 
involving p a vanish thanks to the conservativity of G and we get (|5.57p • Reciprocally, to pass from (15.571) 
to (|2.32p . one has to get back the edge values p a , which can be done exactly as for the hybridization 
of the Generalized Mimetic method in Section [3.11 (using G = G(K,a) such that G(K,(t)k.<j = 1, 
G(K , ct)l,<t = — 1 if L is the control volume on the other side of a and G(K, a)z.e = for other control 
volumes Z and/or edges 8). 

An important element of study of the standard Mimetic method seems to be the existence of a suitable 
lifting operator, re-constructing a flux unknown inside each grid cell from the fluxes unknowns on the 
boundary of the grid cell (see [2] and Section [4~3| . We claim that, for the Generalized Mimetic method 
corresponding to the choice (|5.56[> . the flux F given by [ (|5.52p . (|5.53p ] provides a (nearly) suitable lifting 
operator : it does not completely satisfy the assumptions demanded in [3J Theorem 5.1], but enough so 
that the conclusion of this theorem still holds. 
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Proposition 5.1 For Fk G Fk> let Fk be the restriction to K of F defined by (|5 . 53|) . Then the operator 
Fk G Tk Fk G L 2 (K) satisfies the following properties: 

VFjc e^.Vae^.Viea: Fjf (x) • n A> = F K ,„ , (5.60) 

V-Fr: G Tk , Vg affine function : / g(x)div(Fjf )(x) da; = / q(x)VTV h (FK)wE(x) dx , (5.61) 

Jk Je 

VF G R d , de^mng F K = (F ■ n K , a )ae£ K ■ F K = F , (5.62) 

/or any We satisfying (|2.14|) ■ 

Remark 5.2 Properties (|5.60[) and (|5.62|1 are i/ie same as in |H Theorem 5.1], but Property (|5.61|) zs 
replaced in this reference by the stronger form "div(-FK-) = T)IV h {FK) on K" (xk is also taken as the 
center of gravity, which corresponds to we = I in (|5.61[) ). 

Proof of Proposition [5TT1 

We already noticed (|5.60p (consequence of (|2.30|) and the fact that (x — Xk) ■ t^k,<t — d-K.a for all x G er). 
If F G R d , and F K = {F ■ n K ,a)ae£ K , then ^TM) and (pT2"7"]) show that K k ^k{Fk) = —F and thus that 
T K (F K ) = 0, in which case F K = -A K v K {F K ) = F and (|51)2"]) holds. 
Let us now turn to (|5.6ip . For g = 1, this relation is simply (|5.59[) . If q(x) = x, then 

/ q(x)dxv(F K )(x)dx = V T K ,a(F K )d [ T^dx 



and (|5.54|) then gives 



[ q(x)div(F K )(x)dx = J2 TkAF k )cI^\ A K , 
Jk ±? d K ,* 



cr££ K 



= 2J \ a \ T K,cr{FK)xK 



ae£ K 



^2 W\Fk,ct x k 

\<t<=£k I 



q(x )we(x) dx 



>K 



by assumption on we- ■ 

5.3 Two-point flux cases 

We now consider isotropic diffusion tensors: 



A = \(x)Id, (5.63) 



with X(x) G R and exhibit cases in which the HMMF provides two-point fluxes, in the sense that the 
fluxes satisfy F E = te,e'{pe — Pe 1 ), in the case where e is the common edge of two neighboring cells 
E and E' , with te,e' > only depending on the grid (not on the unknowns). Recall that using a 
two-point flux scheme yields a matrix with positive inverse, and is the easiest way to ensure that, in 
the case of a linear scheme, monotony and local maximum principle hold. The two-point flux scheme is 
also probably the cheapest scheme in terms of implementation and computing cost. Moreover, although 
no theoretical proof is yet known, numerical evidence shows the order 2 convergence [20] of the two- 
point scheme on triangular meshes, taking for xk the intersection of the orthogonal bisectors [T8| [13]. 
Therefore, whenever possible, one should strive to recover this scheme when using admissible meshes (in 
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the sense of [Tj5] previously mentioned). In the HMMF framework, two-point fluxes are obtained, using 
the notations provided by Definition 12.11 if the matrix of the bilinear form [F, G] X h in (|2.6p is diagonal. 
This implies, in the case of meshes such that two neighboring grid cells only have one edge in common, 
that the matrices defining the local inner product (|2.3p are diagonal. If the matrix Me is diagonal, 
we get from the property M^N^ = Re and from (|5.63p that there exists (x e E £ R such that 

fj,%n% =x e — x E . (5.64) 

This implies that xe is, for any face e of E, a point of the orthogonal line to e passing through x e [x E 
is then necessarily unique) and that \f E is the orthogonal distance between xe and e. In the case of 
a triangle, xe is thus the intersection of the orthogonal bisectors of the sides of the triangle, which is 
the center of gravity only if the triangle is equilateral; hence , except in this restricted case, the original 
Mimetic method cannot yield a two-point flux method. Note also that there are meshes such that the 
orthogonal bisectors of the faces do not intersect, but for which nevertheless some "centers" in the cells 
exist and are such that the line joining the centers of two neighboring cells is orthogonal to their common 
face: these meshes are referred to "admissible" meshes in [T»U Definition 9.1 p. 762]; a classical example 
of such admissible meshes are the general Voronoi meshes. On such admissible meshes, the HMMF does 
not provide a two-point flux scheme for isotropic diffusion operators, although a two-point flux Finite 
Volume scheme can be defined, with the desired convergence properties (see [13]). 

In this section, we shall call "super-admissible discretizations" the discretizations which fulfill the property 
(|5.64p for some choice of {xe)e£Si h ■ We wish to show that for all super- admissible discretizations and in 
the isotropic case (|5.63p . the HMMF provides a two-point flux scheme. Using the notations of Definitions 
12.41 and I2.6[ Q5.64j) is written n^.o- = (x ff — xk)/<Ik,<j and defines our choice of parameters (xk )kem- 

Let us take ax,a = Xk in the Hybrid presentation (|2.20p . denoting by Xk the mean value of the function 
X(x) in K. Using Definition (|2.16p for V kPk and thanks to (|2.17p . a simple calculation shows that 

j ^ K S K ,(t(Pk)Sk,AQk) = j ^ K (j*K -Pa){qK - go-) ~ \K\X K V ' kVK ■ Vifgk- 



Hence we deduce from (|2.20p that Fk.o = ~j^{pk ~Pa) and the conservativity (|2.ip leads, for an internal 
edge between 
flux becomes 



edge between the control volumes K and L, to p a — rfA ^ L ^^^'°'^ pg ■ The resulting expression for the 



771 X K Xl , s 

FK ' a = 1 TTT^ y-(PK-PL), 

which is the expression of the flux for the standard 2-points Finite Volume scheme with harmonic averaging 
of the diffusion coefficient. 

It is also easy to find back this expression from the Mixed presentation (|2.3ip . taking 



and, for a ^ a 1 , 



The property 



l K,cr,a 



*K,cr,<T' 



0. 



V \a\^T Kt<T (F K )T K , a (G K ) = V \a\=^-F Kta G K , a -\K\v K (F K )-X K VK(G K ), 

— — *K 

which results from (|2.27p . (I2.30p . (|2.17p and (|5.64|) (under the form d,K,<r^K,a = x a — Xk), shows that 



{Fk,Gk)k — ^ H Fk,(tGk,(j- 
Thanks to (|2.32p . this gives Fx, a = j^(pk — Pa) and we conclude as above. 
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5.4 Elimination of some edge unknowns 

In the study of the hybrid version (|2~25| of the HMMF [T51 [T5], it was suggested to replace the space 
Hm,o by the space Hj^.q defined by H^o = {p ^ Hm,o such that p a = J2kgm„ PvPk if ct £ B}, where: 

i) M. a is a subset of M, including a few cells (in general, less than d + 1, where we recall that d is 
the dimension of the space) "close" from a, 

ii) the coefficients /3jr are barycentric weights of the point x a with respect to the points Xki which 
means that J2i(eM„ $a = 1 and Ekgm„ Pa x K = x„, 

iii) B is any subset of the set of all internal edges (the cases of the empty set or of the full set itself 
being not excluded). 

Then the scheme is defined by: find p £ H%[ such that 

We-ff^o : ]T \ K \^KVKPK-V K qK+ Sk(Qk) T KSk(pk) = ]T q K f f- (5.65) 

KM KeM K<EM K 

In the case where B = 0, then the method belongs to the HMMF and in the case where B is the full 
set of the internal edges, then there is no more edge unknowns, and we get back a cell-centered scheme. 
In the intermediate cases, we get schemes where the unknowns are all the cell unknowns, and the edge 
unknowns p a with a ^ B. 

This technique has been shown in [15] and |16j to fulfill the convergence and error estimates requirements 
in the case of diagonal matrices Bj£. It can be applied, with the same convergence properties, to the 
HMMF with symmetric positive definite matrices Bj^ as in Section [U 



6 Appendix 

6.1 About the Generalized Mimetic definition 

We prove here two results linked with the definition of the Generalized Mimetic method: the existence 
of a weight function satisfying (|2.14| and the equivalence between (|2.15j) and (|2.12| . 

Lemma 6.1 If E is a bounded non-empty open subset of R d and xe £ R d , then there exists an affine 
function we ■ R- d R- satisfying (|2.14|) . 

Proof of Lemma 16.11 

We look for £ £ R d such that we(x) = l + (x — xe) = 1 + (x — xe) t !; satisfies the properties (where xe 
is the center of gravity of E). The first property of (|2.14|) is straightforward since J E (x — xe)c\x = 0, and 
the second property is equivalent to \E\xe + J E x(x — xe) t £,c\x = \E\xe\ since J e xe(x — xe) t dx = 0, 
this boils down to 

(^J (x - x E ){x - x E ) T c\x^j £ = \E\(x E - x E )- (6.66) 

Let Je be the d x d matrix J E (x — x E )(x — xe) T dx: we have, for all 77 £ R d \{0}, JeV ■ 77 = J E ((x — 
xe) ■ n) 2 c\x and the function x — > (x — xe) ■ V vanishes only on an hypcrplane of R d ; this proves that 
JeV • V > for all 77 0. Hence, Je is invertible and there exists (a unique) £ satisfying (|6.66[) . which 
concludes the proof. ■ 

Lemma 6.2 Let [•, -]e be a local inner product on the space of the fluxes unknowns of a grid cell E , and 
let Mb be its matrix. Then [-,-]e satisfies (|2.15[) (with we satisfying (|2.14p ) if and only ifM. E satisfies 
([2~T2"|) (withR E defined by ([2"7Hi|) and (C E ,U E ) defined by ([2~TT)j) and ([2~TT|) ). 
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Proof of Lemma 16.21 

It is known [3] that, for the standard Mimetic method, (|2.4p is equivalent to MeNe = Re with Re 
defined by (|2.9[) and Ne defined in 10[) : similarly, it is quite easy to see that (|2.15p (with we satisfying 
(|2.14|) ) is equivalent to 

M^Nb = Re (6.67) 

with Re defined by (|2.13p : indeed, (|2.15p with q = 1 is simply the definition (|2.2p of the discrete 
divergence operator (because J e we(x) dx = \E\) and, with q(x) = Xj, since J E XjWe(x) dx = \E\(xe)j, 
(|2.15p boils down to 

k E k E 

GIMe^ +^0^^)^ = J^G^\e i \(x ei ) j , 

i=l i=l 

which is precisely G^M^Nb)., = G E '(M.E)j with (Me)] the j-th column of Re- Wc therefore only need 
to compare (f2~T2"|) with (JBTHZJ) ■ 

Let us first assume that satisfies (|2.12p . The generic formula (|2.26p implies 

R T E (N E h = \E\^Eh (6.68) 



and thus 



-^ReA^rINe = R E - (6.69) 



Since C|JNe — by definition of Ce, this shows that Me satisfies (|6.67l) . 

Let us now assume that Me satisfies (|6.67p and let us consider the symmetric matrix 

1 



E - —^e^e^E- (6-70) 



By (|6.67p and (|6.69[) . we have MeNe = 0, and the columns of Ne are therefore in the kernel of Me- 
The definition of Ce shows that the columns of Ne span the kernel of C E : we therefore have ker(C £ ) C 
ker(Ms) and we deduce the existence ofafcEx(fc£ — d) matrix A such that 



-E- 



(6.71) 



By symmetry of Me we have AC E = CeA t and thus AC^Cb = CeA t Ce; but C^Ce is an invertible 
(kE — d) x (kE — d) matrix (since Ce is of rank kE — d) and therefore A = CeA t Ce(CeCe) _1 = CeUe for 
some (kE — d) x (kE — d) matrix Ue. Gathering this result with (|6.70p and (|6.7ip . we have proved that 
satisfies (|2.12p for some Ue, and it remains to prove that this last matrix is symmetric definite positive 
to conclude. By (|2.12p and the symmetry of Me and j^-ReA^R^ we have CeUeCe = CeUeCe and, 
since C E is onto and Ce is one-to-one, we deduce Ue = U^. To prove that Ue is definite positive, we 
use (|2.12p and the fact that Me is definite positive to write 

V£ G kcr(Rl) , £ ? : (V E C E ■ (C T E > 0. (6.72) 

This shows in particular that ker(R^) n ker(C|J) = {0} and thus, since ker(R^) has dimension k E — d, 
that the imag e by Ce of ker(Rg) is ~R kE - d . Equation (IB~72l then proves that Ue is definite positive on 
the whole of R fc « _rf and the proof is complete. ■ 

6.2 An algebraic lemma 

Lemma 6.3 Let X , Y and Z be finite dimension vector spaces and A : X — > Y , B : X — > Z be two linear 
mappings with identical kernel. Then, for all inner product {v}y on Y, there exists an inner product 
{■,-}z on Z such that, for all (x,x') G X 2 , {Bx, Bx'}z = {Ax, Ax'}y ■ 
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Proof of Lemma 16.31 

Let N = kei(A) = ker(£?). The mappings A and B define one-to-one mappings A : X/N — > Y and 
B : X/N — ► Z such that, if x is the class of a;, Ax = Ax and Bi: = Bx. We can therefore work with A 
and -B on X/N rather than with A and B on X, and assume in fact that A and B are one-to-one. 
Then A: X Im(A) and j9 : X — > Im(_B) are isomorphisms and, if {•, >}y is an inner product on we 
can define the inner product {■, -}im(B) 011 Im(_B) the following way: for all z, z' G Im(_B), {z, z'}i m (s) = 
{AB _1 z, AB~ X z'}y (this means that {Bx, Ba/} Im ( B ) = {Ax,Ax'}y for all x, x' G X). This inner 
product is only defined on lm(B), but we extend it to Z by choosing W such that Im(B) © W = Z, 
by taking any inner product {•, -}w on W and by letting {z, z'}z — {zb, ^s}im(B) + { z w, z w}w for all 
z = zb + z' w G Z = Im(£>) © W and z' = z' B + z' w G Z. This extension of {■, -}i m (B) preserves the 
property {Bx, Bx'} z = {Ax, Ax'}y ■ ■ 

Remark 6.4 The proof gives a way to explicitly compute {■, -}z from {•, -}y, A and B: find a supplemen- 
tal space G ofker(A) = ker(B) in X, compute an inverse of B between G and Im(B), deduce {-,-}im(B) 
and extend it to Z by finding a supplemental space of lm(B) . 
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